Effect of Vortex Screening on the Bose Glass to Entangled Liquid Transition of Flux 

Lines in Superconductors 
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We study the phase diagram of flux lines in superconductors with columnar pins. Based on 
numerical exact diagonalisation simulations on small clusters, we get two phases of vortices: A low 
temperature pinned glass with diverging tilt modulus and a high temperature delocalised entangled 
vortex liquid. For random potential disorder we find a new phase transition temperature Tbg from 
a pinned Bose glass to an entangled liquid that reduces with increasing vortex density. This occurs 
primarily because vortices screen the disorder potential and generate an effective weaker random 
potential with increasing vortex density. For a fixed fraction of randomly placed attractive columnar 
pins of strength V p i n , we find a Mott insulating phase when the vortex density exactly matches the 
number of pins (B = -B</>). We also find a transition from a strongly pinned Bose glass for B < B^, 
to a weakly pinned Bose glass for B > B^ as the vortex density is varied. 



PACS Numbers: 74.60. G, 05.30. J 



I. INTRODUCTION 



The interest in equilibrium and dynamical proper- 
ties of flux arrays in type II superconductors originates 
both from their technological importance and the rich 
variety of behavior these systems exhibit under differ- 
ent physical conditions. In a classic paper Abrikosov 
in 1957 showed that for applied fields H such that 
H c i < H < H C 2, the magnetic field penetrates in the 
form of quantized flux tubes (<fo = hc/2e), which in the 
absence of disorder form a triangular lattice. In conven- 
tional low temperature superconductors, this Flux Line 
Lattice (FLL) was believed to exist at all temperatures 
upto H c2 (T). With the discovery of high T c supercon- 
ductors, it was realized that, due to much higher transi- 
tion temperatures, reduced effective dimensionality and 
short coherence lengths, fluctuations play an important 
role in deciding the structure and dynamics of FLLs Q . 
One of the most significant consequences of enhanced 
thermal fluctuations on the FLL is its melting into a flux 
liquid phase via a first order transition ||pj] . 

Disorder arising from vacancies and interstitials, twin 
boundaries, grain boundaries and columnar pins also 
modify the structure and dynamics of the vortex lat- 
tice. The presence of strong static disorder is techno- 
logically relevant, for it leads to effective pinning of vor- 
tices thereby leading to high critical currents ||. It also 
can lead to novel glassy phases such as the vortex glass 
and Bragg glass for the case of random quenched disorder 



Columnar defects i.e. linear damaged tracks in the ma- 
terial caused by heavy ion irradiation have emerged as 
very effective pinning centers [|]|| ■ In order to model the 
effect of columnar pins on the FLL, we exploit the map- 
ping of 3D interacting flux lines onto bosons in (2+l)D 
[jl0|,[llj. In the mapped quantum problem, the columnar 
defects naturally map onto a time-independent random 



potential for bosons Jl2| ]. The irreversibility line can be 
interpreted as a phase transition where the vortex liq- 
uid at high temperatures freezes into a Bose glass (BG) 
of vortices pinned by columnar pins at lower tempera- 
tures. Path integral Monte Carlo simulations [H find a 
low temperature BG with patches of ordered region with 
positional and orientational order which melts into an en- 
tangled defected liquid at high temperatures. Also, when 
the vortex density and defect densities are equal (at the 
matching field B^), each flux line is attached to one pin, 
leading to a Mott insulator (MI) phase . Such a Mott 
insulator has_been observed in magnetization relaxation 
experiments 
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In this paper, we study, using numerical exact diago- 
nalisation on small lattices, the different phases of flux 
lines with columnar pins for various densities of vortices 
and disorder strengths. We consider a lattice of iV-sites 
with N v — H x area/(/)o vortices, interacting via a hard 
core potential. We model disorder in two ways: (a) pin- 
ning disorder where a finite fraction of attractive pins, 
each of strength Vpi n , are placed randomly; and (b) a 
random disorder potential at each point. In case (a), we 
find that an entangled vortex liquid is stable against weak 
pins. For high pinning strengths, a Mott insulator is re- 
alised when the number of vortices is equal to the number 
of pins at the matching field B^. Signatures of a strongly 
pinned Bose glass and a weakly pinned Bose glass are also 
seen as the vortex density is tuned across B^. In case (b), 
we find two phases in the density-disorder plane. At low 
disorder an entangled vortex liquid which localizes into a 
pinned Bose glass with increasing disorder. We find that 
the critical disorder strength required to pin the vortex 
liquid increases with increasing vortex density. This im- 
plies that the temperature required to depin the vortices 
is reduced with increasing fields (see Fig. |t]). 

We organize the paper in the following way. In Sec- 
tion |ll| we give the details of our model. In Section [II A 
we discuss our results for pinning disorder, where we can 
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access different pha ses of vortices including the Mott in- 
sulator. In Section IIIB we discuss our simulations for 



the case where each site has a random disorder potential 
and conjecture an interesting experimental implication of 
our phase diagram. 



II. THE MODEL 

Consider a system of N v flux lines in 3D in a magnetic 
field (B) aligned with the z-axis, described by their 2D 
trajectories Vi{z) as they traverse a sample of thickness 
L with Np columnar pins. Their free energy ]T^ | is given 
by 




1 N v 



(1) 



The first term in Eq. ([!]) is the line tension term with 
tilt modulus ei . The second term denotes the interaction 
energy of all vortex pairs on a constant z-plane, where 
Tij = |r j — Tj\ and V(r) the inter-vortex potential. The 
last term denotes Np columnar pins (|| B), modeled by z- 
independent potential Vp placed on randomly distributed 
positions {p£ m }. 

The classical statistical mechanics of Eq. (|l|) is equiv- 
alent to the quantum mechanics of interacting bosons 
interacting with a potential V(r) in 2D with a random 
static potential Vp(r). The partition function is deter- 
mined by the ground-state energy of a fictitious quantum 
hamiltonian [11 12 Using this mapping, the ther- 
mal fluctuations of the 3D-vortices get mapped onto the 
effective quantum fluctuations of bosons in two spatial 
dimensions and one imaginary time dimension. In this 
mapping, the temperature of the vortex system T plays 
the role of the Planck number h. The bending energy 
ei of the flux lines is equivalent to the mass m of the 
bosons, so that like a lighter particle, a softer flux-line 
can wander more. The length of the vortices is equiva- 
lent to the inverse temperature of bosons so that in order 
to simulate a thick sample (L — > oo) one considers the 
ground state of the quantum hamiltonian given by [fL2| , 

H = -tJ2(4 a 3 + h.c) + Vfa) + E turn (2) 

ij ijtj i 

where, operator Oj (aj) annihilates (creates) a boson at 
site i, m = a\a,i is the number operator and p is the 
chemical potential at site z. i is the measure of quantum 
fluctuations of bosons. //$ has two parts: a uniform part 
(i oc (H — H c i ) which fixes the flux line density and Spi 
which represents site disorder. 

The inter-vortex potential can be written as V(r) = 
2eoK (r/X), with the modified Bessel function K (x) oc 



— \n(x) as x — > 0, and Kq(x) oc a; -1 / 2 exp (— x) for 
x — > oo. The energy scale is set by eo = (0o/4ttA) 2 , 
where 0o = hc/2e and A the London penetration depth. 
In the low density limit i.e. when the inter-vortex sep- 
aration (ao) is large compared A, the long range part of 
V{r) dominates. In the absence of disorder, at low tem- 
peratures, the vortices form a triangular lattice. In the 
high density limit (ao < A), the short ranged part of 
V{r) decides the structure and elastic properties of the 
flux array. In the presence of columnar pins at low den- 
sity, individual vortices get pinned whereas at high den- 
sity, vortex bundles are pinned and the system breaks up 
into patches of ordered regions |14|. We use a simplified 
model and describe the inter boson interaction by a hard 
core interaction. In the low density limit, these bosons 
describe individual vortices whereas in the high density 
limit they describe the renormalized bundles of vortices. 

We fix t to be unity and, disorder 5 fa is introduced 
in two ways: (a) Pinning Disorder : Attractive pins of 
strength V p i n are introduced at randomly chosen 25% of 
the sites, and (b) Box Disorder : Random disorder po- 
tentials at each lattice site chosen either from a uniform 
distribution of width [— A/2, A/2] or from a gaussian dis- 
tribution of width A. We work on a 4 x 4 square lattice 
with periodic boundary conditions, and use exact diago- 
nalization techniques (such as Lanczos fT^]) to compute 
the ground state properties of (||). 

Our aim is to understand the combined effects of the 
strong intervortex hard core interaction and columnar 
disorder on the phases. The bosons described by Eq. (^) 
can exist in three different phases |L3||: (a) Kinetic energy 
dominated superfluid phase for bosons that is equivalent 
to a high temperature entangled flux liquid in the vortex 
problem, where flux lines are delocalised and can hop 
freely from one columnar pin to another, (b) Disorder 
dominated Bose glass phase for bosons with a vanishing 
superfluid density p s , that is equivalent to a Bose glass 
vortex phase with a diverging tilt modulus C44 oc p" 1 . (c) 
Interaction dominated gapped Mott insulator phase for 
bosons corresponds to the phase when the vortex density 
matches the density of columnar pins in the vortex prob- 
lem. As opposed to a BG, the MI is incompressible since 
the vortex density remains locked to the density of pins 
over a finite range of external fields 12 and behaves like 



a Meissner Phase with B = B^ instead of B = 



III. RESULTS 

A. Pinning Disorder 

On a 4 x 4 square lattice, we randomly choose a fixed 
fraction (say 25 %, -/V p i n = 4) of sites to be pins, each of 
strength — V p i n . In Fig. [I], we plot the superfluid density 
p s /p, which is a measure of the stiffness of the phase of 
the order parameter, as a function of density p, for several 
values of V p i n /t. Since we work on a lattice, p s /p — > as 
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p — > 1, and the repulsive local interactions give rise to a 
Mott insulator for any values V p i n . In addition, for large 
VpinA, the superfluid density goes to zero at a special 
Mott insulator point when N v = AT pin = 4; the match- 
ing condition where each columnar pin traps exactly one 
vortex. The distribution P(rii), of local density (n;), fol- 
lows a bimodal distribution, with a peak at and another 
peak at 1. Note that the origin of this Mott insulator (for 
B = B^) is different from the MI obtained in the jam- 
ming limit (p — > 1). This Mott insulator can be viewed 
as a Meissner phase for vortices at the matching field 
instead of zero field P^Jlq] ■ 
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In the SBG phase, vortices are collectively pinned by 
columnar pins 

JUS whereas in the WBG phase, one 
has patches of ordered regions H|. It is interesting to 
see that calculations on such small systems, do exhibit 
signatures of these distinct glassy phases. 
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FIG. 2. Pinning Disorder: Binding energy per vortex 
as a function of density p for N P i„ — 4. Vpin values same 
as in Fig jjj. At large V p i n , vortices are strongly pinned (the 
strongly pinned Bose glass phase) below the matching con- 
dition (p = 0.25) than above (the weakly pinned Bose glass 
phase). 



FIG. 1. Pinning disorder: ps/p as a function of 
p for Vpi n = 0.2,1.0,2.0,3.0,5.0,7.5,8.0 (top to bottom). 
Pins of strength — Vpi n are placed on randomly chosen sites. 
ATpin = 4. Data averaged over 300 disorder realizations. For 
high Vpin the Mott Insulator is realized at the matching field, 
when p — Npi n (indicated by filled triangle on the x-axis) . 

For high Vp- m (say V p - m /t — 8.0 in Fig. |l|), we note that 
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p 



> — 

B>B^ p 



(3) 



implying the vortices are more entangled above the 
matching condition, than below. Alternatively, in Fig. |^ 
we plot binding energy per vortex defined as 



B.E. = 



£(Vp in ) - E(0) 

N v 



(4) 



It is clear that above the matching field (B^) vortices 
are less bound by disorder than below. So depending 
upon the induced B, compared with B^, we find three 
phases within the BG, namely a strongly pinned Bose 
glass (SBG for B < B$), a Mott insulator (for B = 
and a weakly pinned Bose glass (WBG for B > B<j>), 
consistent with ref. Jl9|] . 



B. Box Disorder 



To further illustrate the effect of disorder on the FLL, 
we now present our results for the case where each site i 
on the lattice has a random potential Aj chosen from a 
uniform distribution of [—A/2, A/2]. 

In Fig. U we plot the inverse tilt modulus c^ 1 oc p s as 
a function of A for a fixed density of vortices (p = 0.25). 
For small A, C44 has a finite value (the entangled phase 
of vortices) and for large A it tends to diverge (the Bose 
glass phase of vortices). In the thermodynamic limit, it 
is expected that p s ~ |A — A c \^ with £ = vz in 2D, 
where v is the correlation length exponent, and z is the 
dynamical exponent jl3|]. However, in a finite system, 
the transition is rounded. As an indicator of the critical 
disorder strength A c , we choose the point of maximum 
slope of p s (A)/p (shown on the dotted curve at the bot- 
tom in Fig g). A locus of A c for different densities p 
gives us the phase boundary between the superfluid and 
the Bose glass (see Fig. ||). 

For zero disorder, p s /p differs from unity, which is a 
consequence of broken Galilean invariance. The kinetic 
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energy per particle along the direction of the twist was 
shown to be an upper bound to p s /p |pofl . We plot 
(ke) x /2t/5 as a function of A in the same figure. Near 
A = the kinetic energy per particle is exactly equal 
to the superfiuid density and with increasing disorder it 
serves as an upper bound but being a local quantity, it 
fails to pick up the transition pQ| . 




FIG. 3. Box Disorder: "Superfiuid" density p a j ' p oc cj/ 
(solid triangles) as a function of A for N v — 4; p = 0.25. 
Slope of p 3 (A)/p curve is shown at the bottom (dotted line). 
The point of maximum slope is identified as A c (shown by an 
arrow). The kinetic energy along the direction of the twist, 
(k.e.)x/2tp (solid line), serves as an upper bound to p 3 /p- 
Data averaged over 400 disorder realisations. 

A quantity of interest to probe the spatial extent of the 
wavefunction is the inverse participation ratio defined as 

lHo\Hr)dr 



I = 



(/hM 2 (r)dr) 



(5) 



where, ipo (r) is the ground state wavefunction of Eq. (|^). 
For a localized state, / has a finite value, whereas for an 
extended state, 1^0. In Fig. || we plot 7 as a function of 
A for p = 0.25 and 0.50. 7^0 till A c and starts rising 
above A c indicating localization or pinning of vortices. 
For high disorder, the states are localized for all densi- 
ties. Though it should be noted that, such a numerical 
scheme cannot really distinguish between an extended 
state and a localized state with large localization length 
(say for low A). 

The degree of wandering of a flux line transverse to its 
length due to thermal fluctuations, can be obtained from 
the properties of the ground state wavefunction ipo (r) in 
the mapped problem p2[ by, 



<Tr r 2 %{r) 



(6) 



We plot the transverse wandering length per vortex in 
Fig. | for N v = 5 and 8. Till A c , l± is constant (almost 
equal to the clean case). For A > A c , l± reduces, thereby 
marking localization of vortices due to disorder. 
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FIG. 4. Box Disorder: Inverse Participation Ratio I 
(see Eq.(^J) for definition) which is a measure of the spatial 
extent of the boson wavefunction, for N v = 4 and 8 as a func- 
tion of A. I ~ for entangled vortices (A < A c ) and finite 
for a pinned vortices (A > A c ). 



The locus of critical disorder strength A c , estimated 
from the behaviour of quantities discussed above as a 
function of A, for various densities (p) gives us the phase 
boundary between the superfiuid (entangled liquid) and 
the Bose glass (pinned) in the density-disorder (p — A) 
plane shown in Fig. ^| 

Another measure of the critical disorder strength A c 
(for a given p) is the disorder at which a finite fraction 
of incoherent sites occurs on the lattice. In a superfiuid, 
phase coherence is achieved due to tunneling of bosons 
between different sites leading to fractional occupancies 
on sites. The long wavelength excitations are gapless, 
but due to the phase rigidity, local phase rotations cost 
a finite energy. We call a site coherent if the occupancy 
of that site is fractional. On the other hand, at high 
A, bosons are localised due to disorder. In this disorder 
dominated regime, local phase rotations do not cost any 
energy at a site which is integer occupied. We call a site 
incoherent if (n.;) — or 1, where (n;) = 1 corresponds 
to a single vortex pinned on the site i all along its length. 

A locus of A c (p) based on the emerging of an inco- 
herent site concurs with the phase diagram shown in 
Fig. |^. The phase diagram, is re-entrant due to particle- 
hole symmetry about half filling (p = 0.5). At a mean 
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field level, it has been shown that the system becomes lo- 
calized before clusters of nearly-incoherent sites percolate 
on the lattice Q . 
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FIG. 5. Box Disorder: Transverse wandering length per 
vortex l± (see Eq. (^) for definition) as a function of A for 
N v — 5 and 8. Zj_ is a constant below A c (p), almost equal 
to the clean case, and reduces above A c marking pinning of 
vortices. 

In this mean field theory, it is known that a Bose glass 
can never be accessed, because, a truly incoherent site 
with (rij) = can be obtained only at A = oo pjJ. In 
rcf. ]2Q] , by modeling the disordered boson problem as a 
network of random resistors, it was argued that in order 
to stop the flow of the current (or destroy superfluidity) , 
a line of resistors perpendicular to the direction of the 
flow must be cut. In this perspective, we note that, our 
criterion acts as a precursor to the actual SF-BG tran- 
sition. The true SF-BG transition would lie above our 
A c and below the percolation transition (A£) described 
in ref. @. 

The results presented above are for bounded disorder. 
For the case of unbounded (Gaussian) disorder, we find 
the same qualitative results as above. In some cases, A c 
is marginally reduced relative to the A c for the bounded 
disorder, since one can get large disorder values from the 
tails of the distribution. 

Now we discuss an interesting experimental implica- 
tion of our phase diagram. We see from Fig. || that for 
a higher density of vortices a higher disorder strength is 
required to prevent entanglement and wandering of the 
flux-lines. This happens because with increasing den- 
sity (p), the disorder is effectively screened by vortices 
(strongest attractive sites are occupied first) and a new 
vortex introduced sees a weaker disorder landscape; thus 
is weakly pinned. Therefore, we conjecture that, the tran- 



sition temperature Tbg, for a transition from a pinned 
Bose glass to an entangled liquid, will reduce with in- 
creasing vortex density. 
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FIG. 6. Phase Diagram: Locus of critical disorder A c , 
as a function of vortex density p. The line separates two 
phases: at low disorder, an entangled "superfluid" phase of 
flux lines and at high disorder a pinned "Bose-Glass" phase. 
The separatrix A c (p) is determined by the emergence of a 
finite fraction of incoherent sites with (ni) = or 1 Simu- 
lations are done on a 4 x 4 lattice, with periodic boundary 
conditions. Disorder (Spi) is chosen from a uniform distri- 
bution of width [— A/2, A/2]. Results are averaged over 400 
disorder realizations. 

This conjecture is valid in two regimes in the H — T 
plane, (a) At low fields, where the pinning energy dom- 
inates over interaction energy. For example, when the 
applied field is less than the matching field B^, colum- 
nar pins outnumber vortices and the pinning energy 
(-Spin ~ Co ln(i>/£), where b is the diameter of the colum- 
nar pin and £ the coherence length) dominates over the 
magnetic repulsion fjl| . We show this Tbg by a solid 
line in Fig. ^. The transition line we conjecture has a 
similar shape to the Nelson-Le Doussal line p4|] above 
which point-defects play no role in deciding the struc- 
ture of vortices. The physical reasons for the origin of 
Nelson-Le Doussal line and our Tbg are very different. 
Point defects "promote" flux line wandering and entan- 
glement, whereas correlated disorder "inhibits" wander- 
ing and promotes localization. And it is precisely this 
role of columnar pins which leads to the new Tbg at low 
fields, (b) At high fields, where vortex bundle pinning 
dominates over individual vortex pinning [ p5[ . In our 
model, since we neglect the long-range part in the inter- 
vortex interaction, we do not see signatures of vortex 
bundle pinning at high densities. For the case of vortex- 
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bundle pinning, the screening argument would still hold, 
and hence the topology of Tbg would be the same at high 
fields. 
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FIG. 7. Schematic Phase diagram for vortices in the H — T 
plane. The flux line lattice melts much before H^ F (T). The 
melting line Tbg separates a pinned Bose-Glass and an entan- 
gled vortex liquid. The dashed line is the re-entrant melting 
line. Our conjecture for the Tbg is shown by the solid line, 
in the disorder dominated, low H region. The figure is not 
drawn to scale. 

At intermediate fields, the re-entrant "nose" (shown by 
the dashed line in Fig. ^) was predicted by Nelson jlfl] 
and has been experimentally verified p(J . We dont cap- 
ture this re-entrant part in our simulations since we use 
an approximate form of the inter-vortex potential. 

To conclude, we conjecture a double re-entrant phase 
diagram. Upon decreasing the field near the "nose" the 
vortex liquid localises into a glass where vortex bundles 
are pinned due to disorder and then delocalizes again into 
an entangled liquid. At lower fields, individual vortices 
are pinned by columnar defects, giving rise to another 
disorder dominated glassy phase. 

Finally we note that, the bosonic Hubbard model 
is a paradigm for studying quantum phase transitions 
p7f , in systems such as 4 He on random substrates, su- 
perconducting disordered thin films psj] , magnets and 
Josephson-junction arrays p9| . The phase diagram in 
Fig. ^, for example is consistent with the field induced 
superconductor-insulator (SI) transition observed in cou- 
pled Josephson-junction arrays p9| ]. In low magnetic 
fields, vortices at T = are pinned by disorder intro- 
duced via random offset charges. With increasing fields, 
the vortex density increases and at some critical density, 
vortices Bose-condense |}0| . This vortex superfluid leads 
to an inifinite resistance, giving way to an insulator. 
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